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Abstract 

Graphene is the first example of truly two-dimensional crystals - it's just one layer of carbon 
atoms. It turns out to be a gapless semiconductor with unique electronic properties resulting from 
the fact that charge carriers in graphene demonstrate charge-conjugation symmetry between elec- 
trons and holes and possess an internal degree of freedom similar to "chirality" for ultrarelativistic 
elementary particles. It provides unexpected bridge between condensed matter physics and quan- 
tum electrodynamics (QED). In particular, the relativistic Zitterbewegung leads to the minimum 
conductivity of order of conductance quantum /h in the limit of zero doping; the concept of 
Klein paradox (tunneling of relativistic particles) provides an essential insight into electron prop- 
agation through potential barriers; vacuum polarization around charge impurities is essential for 
understanding of high electron mobility in graphene; index theorem explains anomalous quantum 
Hall effect. 
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Introduction 

The variety of crystallographic forms of carbon places this element into the focus of 
attention both in terms of basic research as well as applications. The tree-dimensional 
crystallographic forms - graphite and diamond - are known from the ancient times, and 
are widely used in industrial applications. Recently discovered zero- dimensional (fullerenes 
or cage moleculesi'^'^) and one- dimensional (carbon nanotubes^) forms are now extensively 
studied due to its' remarkable and, often, unique mechanical and electronic properties. At 
the same time, despite very intensive research in the area, no any two-dimensional form of 
carbon has been known until very recently. 

Ironically, this elusive two-dimensional form (dubbed graphene), is, probably, the best 
theoretically studied carbon allotrope. Graphene - planar, hexagonal arrangements of car- 
bon atoms has been the starting point in all calculations on graphite, carbon nanotubes 
and fullerenes since late 40s^. However, its' experimental discovery has been postponed till 
2004 when a technique called micromechanical cleavage has been employed to obtain first 
graphene crystals^'''. The observation of a peculiar spectrum of charge carriers and anoma- 
lous quantum Hall effect (QHE) in graphene^i^ has initiated enormously growing interest to 
this field (for review, see Refs.— >ii). 

One of the most interesting aspects of the physics of graphene is that it provides a novel 
example of a "feedback" of condensed matter and material science on the fundamental 
physics. Of course, this is not an unique case; archetypical examples are the concept of 
spontaneously broken symmetry playing a crucial role in modern high energy physics and 
quantum field theory^^ and the use of Mossbauer effect to check the general relativity the- 
oryi^. At the same time, such relations are rather rare and always turn out to be very 
fruitful. Actually, discovery of graphene has opened new ways to study some basic quantum 
relativistic phenomena which have always been considered as very exotic. Probably the most 
clear example is the Klein paradox—*^, that is, a property of relativistic quantum particles 
to penetrate with a probability of the order of unity through very high and broad potential 
barriers. Previously it was discussed only for experimentally unattainable (or very hard to 
reach) situations such as particle-antiparticle pair creation at the black hole evaporation, 
or vacuum breakdown at collisions of super-heavy nuclei. At the same time, it appeared 
to be relevant for graphene-based electronics^. Here we discuss similarities and differences 
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FIG. 1: Left: Crystallographic structure of graphene. Atoms from different sublattices (A and B) 
are marked by different shades of gray. Right: Band structure of graphene in the vicinity of the 
Fermi level. Conductance band touches valence band at K and K' points. 

between physics of charge carriers in graphene and quantum electrodynamics (QED). 

I. ELECTRONIC STRUCTURE OF GRAPHENE 

From the point of view of its electronic properties, graphene is a two-dimensional zero-gap 
semiconductor with the energy spectrum shown schematically in Fig. [1] and its low-energy 
quasiparticles formally described by the Dirac-like Hamiltoniani^>^>^ 

Hq = —ikvpcrW (1) 

where vp ~ 10^ ms~^ is the Fermi velocity, and a = [a^, ay) are the Pauli matrices. Neglect- 
ing many-body effects, this description is accurate theoretically^>^>^ and has also been 
proven experimentally^!^ by measuring the energy-dependent cyclotron mass in graphene 
(which yields its linear energy spectrum) and, most clearly, by the observation of a relativis- 
tic analogue of the integer QHE which will be discussed below. 

The fact that charge carriers in graphene are described by the Dirac-like equation ([T]) 
rather than the usual Schrodinger equation can be seen as a consequence of graphene's 
crystal structure, which consists of two equivalent carbon sublattices A and B^>^>^ (see 
Fig. [1]). Quantum mechanical hopping between the sublattices leads to the formation of two 
energy bands, and their intersection near the edges of the Brillouin zone yields the conical 
energy spectrum near the "Dirac" points K and K'. As a result, quasiparticles in graphene 
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exhibit the hnear dispersion relation E = hkvp, as if they were massless relativistic particles, 
with the role of the speed of light played by the Fermi velocity vp ~ c/300. Due to the 
linear spectrum, one can expect that graphene's quasiparticles behave differently from those 
in conventional metals and semiconductors where the energy spectrum can be approximated 
by a parabolic (free-electron-like) dispersion relation. 

Although the linear spectrum is important, it is not the only essential feature that under- 
pins the description of quantum transport in graphene by the Dirac equation. Above zero 
energy, the current carrying states in graphene are, as usual, electron-like and negatively 
charged. At negative energies, if the valence band is not completely filled, its unoccupied 
electronic states behave as positively charged quasiparticles (holes), which are often viewed 
as a condensed-matter equivalent of positrons. Note however that electrons and holes in con- 
densed matter physics are normally described by separate Schrodinger equations, which are 
not in any way connected (as a consequence of the Seitz sum rule^, the equations should 
also involve different effective masses). In contrast, electron and hole states in graphene 
are interconnected, exhibiting properties analogous to the charge-conjugation symmetry in 
QEI>^i^. For the case of graphene, the latter symmetry is a consequence of its crystal 
symmetry because graphene's quasiparticles have to be described by two-component wave- 
functions, which is needed to define relative contributions of sublattices A and B in the 
quasiparticles' make-up. The two-component description for graphene is very similar to the 
one by spinor wavefunctions in QED but the "spin" index for graphene indicates sublattices 
rather than the real spin of electrons and is usually referred to as pseudospin a. 

There are further analogies with QED. The conical spectrum of graphene is the result 
of intersection of the energy bands originating from sublattices A and B (see Fig. [1]) and, 
accordingly, an electron with energy E propagating in the positive direction originates from 
the same branch of the electronic spectrum as the hole with energy —E propagating in 
the opposite direction. This yields that electrons and holes belonging to the same branch 
have pseudospin a pointing in the same direction, which is parallel to the momentum for 
electrons and antiparallel for holes. This allows one to introduce chirality^, that is formally 
a projection of pseudospin on the direction of motion, which is positive and negative for 
electrons and holes, respectively. The term "chirality" is often used to refer to the additional 
built-in symmetry between electron and hole parts of graphene's spectrum and is analogous 
(although not completely identicaP^'^) to the chirality in three-dimensional QED. 
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An alternative view on the origin of the chirahty in graphene is based on the concept 
of "Berry phase" Since the electron wave function is a two-component spinor, it has to 
change sign when the electron moves along the close contour. Thus the wave function gains 
an additional phase (p = tt. 

The analogy with the field theory takes a very interesting twist, should we take into 
account that a sheet of graphene must always be corrugated. The fact that graphene crystals 
always exhibit some finite local curvature can be considered as a consequence of the Mermin- 
Wagner theorem, and has been confirmed experimentally both for graphene samples resting 
on a substrate^^, as well as for free-hanging graphene films^^. 

It is well known that harmonic approximation in the two-dimensional case doesn't pro- 
duce a solution with long-range order ^i2^i21i2^. One can see this as bending instabilities, 
due to soft long-wavelength phonons, lead to membrane crumpling^. Anharmonic coupling 
between bending and stretching modes changes the situation drastically and prevents the 
crumpling^i^'^. However, the membrane should be rippled in a sense that typical fluctua- 
tions in the direction perpendicular to the surface h {x, y) has a scale of order of a (L/a)'' ^ a 
where a is the lattice constant, L is the size of the sample and C, is the roughness exponent. 
The latter can be estimated as C ^ C.G^i^. For a typical sample size L ~ l/xm, the typical 
amplitude of the corrugation for free-hanging membrane at room temperature was estimated 
as 0.5 nm, with their characteristic size being about 5nm— . 

These ripples lead to important consequences for the electronic structure of graphene. 
The nearest-neighbor hopping integral 7 turns out to be fluctuating due to its dependence 
on the deformation tensor— 

_ 1 / dui ^ duj ^ duk duk ^ dh dh\ . . 

2\ dxj dxi dxi dxj dxi dxj J 

where Xj = (x, y) are coordinates in the plane and Ui are corresponding components of the 
displacement vector: 

Taking into account this inhomogeneity in a standard tight-binding description of the elec- 
tronic structure of graphene^ one can obtain, instead of Eq.([T]) an effective Dirac-like Hamil- 
tonian describing electron states near the i^-point: 

H = VF<y {-ihV - -^A^ (4) 

5 



where vp = V3'yoa/2h and A is the "vector potential" connected with the deviations of the 
hopping parameters 7^ from their unperturbed value 70: 

■Ax = 7^ — (72 + 73-271), 
zevp 

A -^(73-7.). (5) 

where the nearest neighbors with vectors (—a/ -\/3, O) ; (a/2-\/3, —a/2) ; (a/2-\/3, a/2) are la- 
belled 1,2, and 3, correspondingly. This means that the roughness fluctuations acts on the 
electronic structure near the /T-point as an Abelian gauge fieldr^ which is equivalent to the 
action of random magnetic field. This means that the bending of graphene violates the 
time-reversal symmetry for a given valley; of course, the Umklapp processes between K and 
K' points will restore this symmetry. As was suggested in Ref.— these effective magnetic 
fields are responsible for suppression of the weak localization effects in graphene. 

Whereas smooth deformation of the graphene sheets produces gauge field similar to elec- 
tromagnetic one, different topological defects in graphene inducing inter-valley (Umklapp) 
processes can be considered as sources of a non- Abelian gauge field; corresponding analogy 
with gravitation was discussed in Refs.— 

In relativistic quantum mechanics, chirality is a consequence of particle-antiparticle sym- 
metry, which also guarantees linear energy spectrum for massless particles. The discovery 
of graphene opens a unique opportunity, to investigate chiral particles with parabolic ( non- 
relativistic) energy spectrum. Quasiparticle with such unusual properties can be found in 
bilayer graphene^. For two graphene layers, the nearest-neighbor tight-binding approxi- 
mation predicts a gapless state with parabolic touching in K and K' pointa^"^ (instead 
of conical band crossing in graphene). The electronic spectrum in this approximation is 
described by a single-particle Hamiltonian^^i^ 

j -{px-ipyf /2m\ 

y-iPx + iPyf /^m J 

where pi = —ihd/dxi are electron momenta operators and m is the effective mass. Here 
we neglected higher-order hopping processes which are important only for very low Fermi 
energies. More accurate consideration:^ gives a very small band overlap (about 1.6 meV) 
but at larger energies bilayer graphene can be treated as a gapless semiconductor. At the 
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same time, electronic states are still characterized by chirality and by non-trivial Berry phase 
2yj-35^ (^in contrast to the case of graphene, where the Berry phase was found to be tt—^). 

II. CHIRAL TUNNELING AND THE KLEIN PARADOX 

The term Klein paradox— li^ii^ usually refers to a counterintuitive relativistic process in 
which an incoming electron starts penetrating through a potential barrier if its height Vq 
exceeds twice the electron's rest energy mc^ (where m is the electron mass and c the speed 
of light). In this case, the transmission probability T depends only weakly on the barrier 
height, approaching the perfect transparency for very high barriers, in stark contrast to 
the conventional, nonrelativistic tunneling where T exponentially decays with increasing 
Vq. This relativistic effect can be attributed to the fact that a sufficiently strong potential, 
being repulsive for electrons, is attractive for positrons and results in positron states inside 
the barrier, which align in energy with the electron continuum outside. Matching between 
electron and positron wavefunctions across the barrier leads to the high-probability tunneling 
described by the Klein paradox. 

One can think about non-relativistic quantum mechanical tunneling though a potential 
barrier in terms of indeterminacy principle. Since momentum and coordinates of a particle 
can not be measured simultaneously a particle can propagate through a classically forbidden 
region where the momentum is formally imaginary, and only coordinates are well defined. 
In relativistic quantum mechanics even coordinate itself cannot be measured with arbitrary 
accuracy, due to pair creation at this measurement. In another words, the Klein paradox 
illustrates that the relativistic quantum mechanics can be consistently formulated only in 
terms of fields rather than individual particles^^. 

Although Klein's gedanken experiment is now well understood, the notion of paradox is 
still used widely, perhaps because the effect has never been observed experimentally. Indeed, 
its observation requires a potential drop ~ mc^ over the Compton length h/mc, which yields 
enormous electric fields^'^(£^ > lO^^V/ cm) and makes the effect relevant only for such exotic 
situations as, for example, positron production around super-heavy nuclei^'^ with charge 
Z > 170 or evaporation of black holes through generation of particle-antiparticle pairs near 
the event horizon^^. At the same time, electronic structure of graphene provides us an 
unique opportunity for easy experimental realization of the Klein ultrarelativistic tunneling 
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in p — n junctionsi^i^i^. 

Let us consider for simplicity the potential barrier that has a rectangular shape and is 
infinite along the y-axis: 

( Vo, 0<x <D, 
V{x)={ (7) 
I otherwise. 

This local potential barrier inverts charge carriers underneath it, creating holes playing the 
role of positrons. For simplicity, we assume in Eq.([7]) infinitely sharp edges, which allows 
a direct link to the case usually considered in QED^^. The sharp-edge assumption is 
justified if the Fermi wavelength A of quasiparticles is much larger than the characteristic 
width of the edge smearing, which in turn should be larger than the lattice constant (to 
disallow Umklapp scattering between different valleys in graphene). Such a barrier can be 
created by the electric field effect using a thin insulator or by local chemical doping^i^i^. 
Importantly, Dirac fermions in graphene are massless and, therefore, there is no formal 
theoretical requirement for the minimal electric field S to form positron-like states under 
the barrier. To create a well-defined barrier in realistic graphene samples with a disorder, 
fields S ~ lO^V/ cm routinely used in experiments-'- should be sufficient, which is eleven 
orders of magnitude lower than the fields necessary for the observation of the Klein paradox 
for elementary particles. 

It is straightforward to solve the tunneling problem for Dirac electrons^^. We assume 
that the incident electron an angle with respect to the and then 

try the components of the Dirac spinor ipi and ip2 for the Hamiltonian H = Hq + V {x) in 
the following form: 



^1 {x,y) 



^2 {x,y) 



(ae*"?^^ + 6e-''?^^) e'^^y , < x < D, 

s' (ae*^^^+^^ - 6e-^«-^-^^) e^'^^y, < x < D, (8) 

g-j^^ikxX+ikyy+icj) X > D 



where kp = 27t / X is the Fermi wavevector, k^ = kp cos and ky = kp sin (p are the wavevec- 
tor components outside the barrier, = \J {E — Vof /h?Vp — fc^, 6 = tan~^ (ky/qx) is the 
refraction angle, s = signE, s' = sign{E — Vq). Requiring the continuity of the wave- 
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FIG. 2: Transmission probability T through a 100-nm-wide barrier as a function of the incident 
angle for single- (a) and bi-layer (b) graphene. The electron concentration n outside the barrier is 
chosen as 0.5 x 10 cm^ for all cases. Inside the barrier, hole concentrations p are 1 X 10^2 and 
3 X 10^^ cm~^ for solid and dashed curves, respectively (such concentrations are most typical in 
experiments with graphene). This corresponds to the Fermi energy E of incident electrons ~ 80 
and 17 meV for single- and bi-layer graphene, respectively, and A ~ 50 nm. The barrier heights 
are (a) 200 and (b) 50 meV (solid curves) and (a) 285 and (b) 100 meV (dashed curves). 

function by matching up coefficients a, 6, t, r, we can find the the reflection coefficient r. 

Fig. [2^ shows examples of the angular dependence of transmission probability T = \t\^ = 
1 — |r|^ calculated using the above expression. The barrier remains always perfectly trans- 
parent for angles close to the normal incidence = 0. The latter is the feature unique to 
massless Dirac fermions and directly related to the Klein paradox in QED. One can un- 
derstand this perfect tunneling in terms of the conservation of pseudospin. Indeed, in the 
absence of pseudospin-fiip processes (such processes are rare as they require a short-range 
potential, which would act differently on A and B sites of the graphene lattice), an electron 
moving to the right can be scattered only to a right-moving electron state or left-moving 
hole state. The matching between directions of pseudospin a for quasiparticles inside and 
outside the barrier results in perfect tunneling. In the strictly one-dimensional case, such 
perfect transmission of Dirac fermions has been discussed in the context of electron trans- 
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port in carbon nanotubes^i^. This is also the case for the two-dimensional (2D) problem of 
graphene which can be directly proven by consideration of the Lippmann-Schwinger equation 
for scattering of arbitrary scalar potential (see Ref.— , supplementary information). 

To elucidate which features of the anomalous tunneling in graphene are related to the 
linear dispersion and which to the pseudospin and chirality of the Dirac spectrum, it is 
instructive to consider the same problem for bilayer graphene^. Although "massive chiral 
fermions" with parabolic spectrum near the touching point and Berry phase 27r do not exist 
in the field theory their existence in the condensed matter physics (confirmed experimen- 
tally^^) offers a unique opportunity to clarify the importance of chirality in the relativistic 
tunneling problem described by the Klein paradox. 

Charge carriers in bilayer graphene are described by an off-diagonal Hamiltonian 
which yields a gapless semiconductor with chiral electrons and holes having a finite mass 
m. An important formal difference between the tunneling problems for single- and bi- 
layer graphene is that in the latter case there are four possible solutions for a given energy 
E = ±h?kp/2'm. Two of them correspond to propagating waves and the other two to 
evanescent ones. Accordingly, for constant potential Vi, eigenstates of Hamiltonian (Ej) 
should be written as 

ipi (x, y) = (oie*''"^ + fejC"'''"^ + Qe"'"'' + die'''"''') e*''^^ 

ip2 (x, y) = Si (^aic''"'^^^^''''^ + 6.e-*'=-^-2^<^» - q/i.c'^'-^ - ^e"''"^'^ e^''^^ (9) 

where 

Si = sign {Vi, — E) ; hki^ = \/lm \E — Vi\ cos</)j; hkiy = ^/2/m\W^-\^ sin (pi 



i^ix = \ kf^ + : hi = { \ 1 + sin (hi - sin 



To find the transmission coefficient through barrier ([7]), one should set rfi = for x < 0, 
63 = C3 = for X > D and satisfy the continuity conditions for both components of 
the wavefunction and their derivatives. For the case of an electron beam that is incident 
normally (0 = 0) and low barriers Vq < E (over-barrier transmission), we obtain ipi = —'ip2 
both outside and inside the barrier, and the chirality of fermions in bilayer graphene does 
not manifest itself. In this case, scattering at the barrier ((Tj) is the same as for electrons 
described by the Schrodinger equation. However, for any finite (even in the case Vq < E), 
waves localized at the barrier interfaces are essential to satisfy the boundary conditions. 
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The most intriguing behavior is found for Vq > E, where electrons outside the barrier 
transform into holes inside it, or vice versa. Examples of the angular dependence of T in 
bilayer graphene are plotted in Fig. ^jp. They show a dramatic difference as compared with 
the case of massless Dirac fermions. There are again pronounced transmission resonances at 
some incident angles, where T approaches unity. However, instead of the perfect transmission 
found for normally-incident Dirac fermions (see Fig. [2^), our numerical analysis has yielded 
the opposite effect: Massive chiral fermions are always perfectly reflected for angles close to 
= 0. At the same time, there is always a "magic angle" when the transmission probability 
is equal to one. 

The fact that a barrier (or even a single p — n junction) incorporated in a bilayer graphene 
device should lead to exponentially small tunneling current can be exploited in developing 
graphene-based field effect transistors (FET)^^. Such transistors are particularly tempting 
because of their high mobility and ballistic transport at submicron distances^'^'^. However, 
the fundamental problem along this route is that the conducting channel in single-layer 
graphene cannot be pinched off due to the Klein paradox (alternative view on this fact 
based on the concept of the minimal conductivity will be discussed in the next section). 
This severely limits achievable on-ofF ratios for such FETs^ and, therefore, the scope for 
their applications. A bilayer FET with a local gate inverting the sign of charge carriers 
should yield much higher on-off ratios. 



III. PROBLEM OF MINIMUM CONDUCTIVITY 

One of amazing properties of graphene is its finite minimal conductivity which is of the 
order of the conductance quantum e^/h per valley per spin; it is important to stress that 
this is the "quantization" of conductivity rather than of conductance^. This is not only very 
interesting conceptually but also important in the view of potential applications of graphene 
for ballistic field-effect transistora^*^. At the same time, this phenomenon is intimately 
related with specifically quantum-relativistic phenomenon known as Zitterbewegung^^i^. 

Numerous considerations of the conductivity of a two-dimensional massless Dirac fermion 
gas do give this value of the minimal conductivity with accuracy of some factor of the order of 
unity^'^'^'^'^'^'^'^'^i^. It is really surprising that in this case there is a final conductivity 
for an ideal crystal, that is, without any scattering processes^'^'^'^. 
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The Dirac Hamiltonian ([T]) in the secondary-quantization representation takes the form 

H = vJ2K^P'^r> (10) 
p 

and the corresponding expression for the current operator 

j = et; ^ ^J,(T^p (11) 



where p is the momentum and = ^V'pi)V'p2) pseudospinor electron operators. Here 
we omit spin and valley indices (so, keeping in mind applications to graphene, the results 
for the conductivity should be multiplied by 4 due to two spin projections and two valleys). 
Straightforward calculations result in the following time evolution of the electron annihilation 
operator 

^P it) = 



P 



1 - 



pa 



P 



(12) 



and for the current operator 



j(t) = jo(t)+ji(t)+jI(t) 
t P (P^) . 



jo it) 
31 it) 



p 

2 ^ P 



P 



a 



2 *P 



p2 



+ -a X p 

P 



(13) 



where ep = vp/h is the particle frequency. The last term in Eq. (fT3|) corresponds to the "Zit- 
terbewegung" , a phenomenon connected with the uncertainty of the position of relativistic 
quantum particles due to the inevitable creation of particle-antiparticle pairs at the position 
measurement^. 

To calculate the conductivity a (uj) following Ref.— we will try first to use the Kubo 
formulai^ which reads for two-dimensional isotropic case: 

oo /? 

^ (^) = ^ 1 dte'-' I dX (j {t - a) j) (14) 



where (3 = is the inverse temperature, A is the sample area. In the static limit uj = 
taking into account the Onsager relations and analyticity of the correlators (j(-2)j) for 
—(3 < Imz < one has^ 

oo 
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Usually, for ideal crystals, the current operator commutes with the Hamiltonian and thus j (t) 
does not depend on time. In that case, due to Eq. ffT^ the frequency-dependent conductivity 
contains only the Drude peak 



Either the spectral weight of the Drude peak is finite and, thus, the static conductivity is 
infinite, or it is equal to zero. It is easy to check that for the system under consideration the 
spectral weight of the Drude peak is proportional to the modulus of the chemical potential 
(cf. Eq.(44) of Ref.— ) and thus vanishes at zero doping (yU = 0). It is the Zitterbewegung, 
i.e. the oscillating term ji (t) which is responsible for nontrivial behavior of the conductivity 
for zero temperature and zero chemical potential (that is, the case of no charge carriers). A 
straightforward calculation gives a formal result 



where one delta-function originates from the integration over t in Eq.( fT5ll and the second 
one - from the derivative of the Fermi distribution function appearing at the calculation of 
the average over product of Fermi-operators. Of course, the square of the delta function 
is not a well-defined object and thus Eq. ffTTl) is meaningless before specification of the way 
how one should regularize the delta-functions. After regularization the integral in Eq. (fT7|) 
is finite, but its value depends on the regularization procedure (see Refs.— i^). 

To circumvent the problem of ambiguity in the expression for cr in Eq. lfTTll it is instructive 
to follow the alternative Landauer approach to calculate the conductivity^"^. Let us assume 
that our sample is the ring of length Ly in y direction; we will use the Landauer formula to 
calculate the conductance in the x direction (see Fig. [3]). There is still an uncertainty in 
the sense that the conductivity turns out to be dependent on the shape of the sample. To 
have a final transparency we should keep finite. On the other hand, periodic boundary 
conditions in y direction are nonphysical and we have to choose Ly as large as possible to 
weaken their effects. Thus, for two-dimensional situation one should choose Lx Ly. 

In the coordinate representation the Dirac equation at zero energy takes the form 




(16) 



oo 




(17) 







{Kx + iKy)ilji = 
{Kx - iKy)^lj2 = 



(18) 
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FIG. 3: Geometry of the sample. The thick arrow shows the direction of current, ipt (sohd Hne) 
and ipb (dashed hne) are wave functions of the edge states locahzed near the top and the bottom 
of the sample, correspondingly. 

where Ki = —i-£r-- General solutions of these equations are just arbitrary analytical (or 
complex conjugated analytical) functions: 



ipi = i!i{x + iy) , 
i'2 = V'2 (a; - iy) . 



(19) 



Due to periodicity in the y direction both wave functions should be proportional to exp {ikyy) 
where ky = 2'n'n/Ly, n = 0, ±1, ±2, .... This means that the dependence on x is also fixed: the 
wave functions are proportional to exp {±2iTnx/Ly) . They correspond to the states localized 
near the bottom and the top of the sample (see Fig. [3]). 

To use the Landauer formula, we should introduce boundary conditions at the sample 
edges (x = and x = L^). To be specific, let us assume that the leads are made of doped 
graphene with the potential Vq < and the Fermi energy Ep = vkp = — Vq. The wave 
functions in the leads are supposed to have the same y-dependence, that is, iIji^2 {x, y) = 
i'1,2 (x) exp (ikyy). Thus, one can try the solution of the Dirac equation in the following 
form: 



tpi (x) 



1 



te 



X < 

< X < La, 
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^2 (x) 



be~''y'', 0<x <L^ (20) 



where sm0 = ky/kp,kj. = \/kp — k'l. /^From the conditions of the continuity of the wave 



y 

functions, one can find the transmission coefficient 

Tn=\t{k,)\'= ,,,;f.'^ . , ■ (21) 

cosh [kyLx) — sm (p 

Furthermore, one should assume that kpL^ ^ 1 and put ~ in Eq. fl2T]) . Thus, the trace 
of the transparency which is just the conductance (in units of e^/h) is 

oo 

TrT = y 2^ ~ (22) 

^f-^cosh2(A;,L,.) vrL, 

Assuming that the conductance is equal to one finds the contribution to the conductivity 
equal to e^/ij^h). Experimentally^, it is close to e^//i, that is, roughly, three times larger 
than this estimation. This discrepancy known as a "missed pi(e) problem" is still not solved. 

One should stress that this is the gapless character of the spectrum and the chirality 
rather than the linear dispersion relation which lead to the minimal conductivity. Different 
estimations for the bilayer graphene give the same order of magnitude for the conductiv- 
ityi^iSEi^i^i^i^^ in agreement with experimenatl observations^. 

Another approach to qualitative understanding of the minimal conductivity is based 
on the Klein paradox—. In conventional 2D systems, strong enough disorder results in 
electronic states that are separated by barriers with exponentially small transparency^"^. 
This is known to lead to the Anderson localization. In contrast, in both graphene materials 
all potential barriers are relatively transparent (T ^ 1 at least for some angles) which does 
not allow charge carriers to be confined by potential barriers that are smooth on an atomic 
scale. Therefore, different electron and hole "puddles" induced by disorder are not isolated 
but effectively percolate, thereby suppressing localization. In the absence of the localization, 
a standard Mott estimation of the minimal conductivity assuming that the mean-free path 
cannot be smaller that electron wave-length^ immediately gives the minimal conductivity 
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IV. INDEX THEOREM AND ANOMALOUS QUANTUM HALL EFFECT 



The anomalous QHE in graphene^i^ is probably the most striking demonstration of the 
massless character of the charge carrier spectrum in graphene. In two-dimensional systems 
with a constant magnetic field B perpendicular to the system plane the energy spectrum is 
discrete (Landau quantization). For the case of massless Dirac fermions the energy spectrum 
takes the form^i^ 

= ^J2\e\Bhv^p{u + l/2±l/2) (23) 

where vp is the electron velocity, u = 0, 1,2... is the quantum number and the term with 
±1/2 is connected to the chirality. Just to remind that in the usual case of the parabolic 
dispersion relation the Landau level sequence is E = hvciv + 1/2) where Uc is the frequency 
of electron rotation in magnetic field (cyclotron frequency )^>^. 

An important peculiarity of Landau levels for massless Dirac fermions is the existence of 
zero-energy states (with z/ = and the minus sign in equation (!23|) ). This situation differs 
fundamentally from usual semiconductors with parabolic bands where the first Landau level 
is shifted by tujJc/2. The existence of the zero-energy Landau level leads to an anomalous 
QHE with half-integer quantization of the Hall conductivity, instead of the integer one (for 
a review of the QHE, see, e.g., Ref.~). Usually, all Landau levels have the same degeneracy 
(a number of electron states with a given energy) which is just proportional to a magnetic 
flux through the system. As a result, plateaus in the Hall conductivity corresponding to the 
filling of the first v levels are just integer (in units of the conductance quant e'^/h). For the 
case of massless Dirac electrons, the zero-energy Landau level has twice smaller degeneracy 
than any other level (it corresponds to the minus sign in the equation (1) whereas the energy 
level proportional to with integer p > 1 are obtained twice, for z/ = p and minus sign, 
and for u = p — 1 and plus sign). A discovery of this "half- integer QHE"-i^ was the most 
direct evidence of the Dirac fermions in graphene. 

A deeper view on the origin of the zero-energy Landau level and thus the anomalous 
QHE is provided by the famous Atiyah-Singer index theorem which plays an important role 
in the modern quantum field theory and theory of superstrings^^. The Dirac equation has a 
charge-conjugation symmetry between electrons and holes. This means that for any electron 
state with a positive energy E a corresponding conjugated hole state with the energy —E 
should exist. However, the states with zero energy can be, in general, anomalous. For 
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curved space (e.g., for a deformed graphene sheet with some defects of crystal structure) 
and/or in the presence of so called "gauge fields" (the usual electromagnetic field provides 
just the simplest example of these fields) sometimes an existence of states with zero energy 
is guaranteed by topological reasons, this states being chiral (in our case this means that 
depending on the sign of the magnetic field there is only sublattice A or sublattice B states 
which contribute to the zero-energy Landau level). This means, in particular, that the 
number of these states expressed in terms of total magnetic flux is a topological invariant 
and remains the same even if the magnetic field is inhomogeneous^. This is an important 
conclusion since, as discussed in section 1, the ripples on graphene create strong effective 
inhomogeneous magnetic fields leading, in particular, to suppression of weak localization^^. 
However, due to these topological arguments they cannot destroy the anomalous QHE in 
graphene. About applications of the index theorem to two-dimensional systems, and, in 
particular, to graphene see also Refs.— i^. 

An alternative view on the origin of the anomalous QHE in graphene is based on the con- 
cept of the "Berry phase" which was discussed above in connection with the Klein paradox. 
When chiral electron moves along the close contour its wave function gains an additional 
phase = TT. In quasiclassical terms (see, e.g., Refs.— i^), stationary states are nothing but 
electron standing waves and they can exist if the electron orbit contains, at least, half of the 
wavelength. Due to the additional phase shift by the Berry phase, this condition is satisfied 
already for the zeroth length of the orbit, that is, for zero energy. Other aspects of the QHE 
in graphene are considered in papera^^'^'^'^. 

In the previous section we have discussed an absence of localization in graphene as one 
of the main consequences of relativistic quantum effects such as the Klein paradox. At the 
same time, the localization is of crucial importance for the QHE providing a well-defined 
Hall plateau^. Here we demonstrate that the localization by a scalar potential V in the 
magnetic field is possible. For simplicity, let us assume the case of a weak one-dimensional 
inhomogeneity with the potential V = V{y) much smaller than the cyclotron quantum of 

1 /2 

the order of hvp/l where / = {hc/\e\B) ' is the magnetic length. 

Using a standard gauge for the vector potential = = 0, Ay = Bx one can write the 
Dirac Hamiltonian for the problem as 

H = Hi + H2 (24) 



17 



where 

-^f^^F (I: + f ) 

( Viy) hvF^\ , , 

i^2= ^ \ . (26) 

Let us try eigenfunctions of the Hamiltonian (I24p with the energy E as an expansion in the 
basis of appropriate harmonic oscillator problem: 



u=0 



Vu{ky,x) = (^^^^^^-J^ ] , (27) 

where i = 1,2 is the pseudospinor index and Di,{z) ~ exp {—z^ /A) (^z/\/2j the We- 
ber functions^-. After straightforward calculations one obtains a set of equations for the 
expansion coefficients Cu^ (ky): 

- ^ (1 - ^.o) c(^) (A:,) = (ky) - ^ E / ^^fc. - %) (qy) {i^, ky, qy) , 



-(l + .)c»'(*») = ^ 



z/, fcj^lz/', gy) = / dxify {ky, x) Lp^i {qy, x) , (28) 



v{q) is a Fourier component of l^(y). For the case of a weak (y)| ^ hvp/l) and smooth 
potential one can neglect a mixing of different Landau bands, similar to the case of QHE 
for conventional electron gas^^. As a result of direct calculations one can demonstrate that 
in this case the solution of Eq. fl28l) describes the Landau wave functions of nonperturbed 
problem with the orbit center satisfying the equation 

E±^V2^=V{y,). (29) 

Thus, the Landau levels in the case under consideration are just smeared into the bands of 
localized states with the localization radius of order of the magnetic length, similar to the 
case of the usual QHE^. 
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V. NONLINEAR SCREENING OF CHARGE IMPURITIES 



In QED, vacuum polarization effects due to creation of virtual electron-positron pairs 
by an external potential are of a crucial importances*^. Due to the smallness of the fine 
structure constant oqed = e^//ic ~ 1/137 these effects leading to logarithmic corrections 
to the Coulomb potential are rather small, except the case of very small distances. In 
graphene, the effective fine structure constant a = e'^/hvpe is of the order of unity (e is the 
dielectric constant due to substrate; for the case of graphene on quartz one should choose^ 
e ~ 2.4 — 2.5) and the corresponding effects can be of crucial importance for transport 
properties^i. Due to the absence of explicitely small parameters describing correlation effects 
in graphene it is very difficult to consider this problem rigorously and the situation still is 
controversial. Nevertheless, the vacuum polarization effects are in our opinion important 
and should be taken into account in any future theory. 

A general nonlinear theory of screening in the system of interacting particles can be 
formulated in a framework of the density functional approacb^^. In this theory a total 
potential V (r) acting on electrons is equal 



where Vq (r) is an external potential and Vind (r) is a potential induced by a redistribution 
of electron density: 



where the first term is the Hartree potential and the second one is the exchange-correlation 
potential. We consider here only a redistribution of charge carriers in the external impurity 
potential 



taking into account contributions of the crystal lattice potential and of electrons in com- 
pletely filled bands via a dielectric constant e and compensated homogeneous charge density 
—en. Here Z is the dimensionless impurity charge (to be specific, we will assume Z > 0; it 
can be easily demonstrated that, actually, in our final expressions Z should be just replaced 
by \Z\). This kind of approach is valid at a spacial scale much larger than a lattice constant; 
in all other aspects, it is formally exact until we specify the expressions for V^c and n [V (r)]. 




(30) 




(31) 




(32) 



er 
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The Thomas-Fermi theory^ which is, actually, the simplest approximation in the density 
functional approach, was used in Ref.*"^. It is based on the two assumptions: (i) we neglect 
the exchange-correlation potential in comparison with the Hartree potential in Eq. fl3ip and 
(ii) we put n (r') = — V (r')], where n (/i) is a density of a homogeneous electron gas 
with chemical potential /x. The former assumption means that we are interested in the 
long-wavelength response of the electron system and thus the long-range Coulomb forces 
dominate over the local exchange-correlation effects. The latter one holds provided that the 
external potential is smooth enough. A rigorous statement is that an addition of a constant 
potential is equivalent to the shift of the chemical potential. In particular, the Thomas- Fermi 
theory gives an exact expression for a static inhomogeneous dielectric function e(g) in the 
limit of small wavevectors q —>■ 0^. The Thomas-Fermi theory of atoms is asymptotically 
exact in the limit of infinite nuclear charge^. 

In the Thomas- Fermi theory Eq. (!3T]) reads 

V'^.W^H^/dr'^ii^i-^M^. (33) 

e J \r — r I 

The function n (/i) is expressed via the density of states (E): 

A* 

n{fi) = j dEf (E) N{E) = j dEN (E) (34) 

where / {E) is the Fermi function, and the last equality is valid for zero temperature (we 
will restrict ourselves here only to this case). For the case of graphene with the linear energy 
spectrum near the crossing points K and K' one has 

"("'^^^ (-) 

where we have taken into account a factor 4 due to two valleys and two spin projections. 

Let us start first with the case of zero doping (/i = 0) where, according to the linear 
response theory, there is no screening at all. Substituting Eqs.( l33|) . (!32l) . and (!35|) into 
Eq. (l30l) . introducing the notation 

V (r) = —E (r) (36) 

er 

and integrating over the polar angle of vector r', we obtain a nonlinear integral equation for 
the function F (r) : 



oo 



F(r)^Z-?^ f'^^K ?^]FHr') (37) 
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where 

7r/2 

Kik)= , . , (38) 

J V 1 — fc"' sm ip 



is the complete elhptic integral, 



for the case of graphene on Si02 Q ~ 2. 

An approximate solution of the integral equation (!37|) at r ^ a gives the following result^ 

F (r) ~ (40) 

^ ^ 1 + In ^ ^ ^ 

which is formally similar to that in QED^ but with a much larger value of the interaction 
constant. Note that logarithmically divergent corrections to the electron effective mass in 
undoped graphene due to the vacuum polarization effects have been considered in Ref.— . 

For the case of doped graphene the main result of Ref.— is the replacement of the bare 
charge impurity Z by 

Z* ^ ^. (41) 

1 + ZQln-i ^ ^ 



Ka 

where 

.2 ,,,, 

(42) 



4e^ l/i 

K - 



is the inverse screening radius. This weakens essentially this scattering mechanism since 
Q In ^ is of order of ten for typical charge carrier concentrations. Perturbative estimations 
of the electron mobility®- should be thus multiplied by this factor squared. As a result, the 
mobility for the same parameters turns out to be two order of magnitude larger. Instead of 
a concentration-independent mobility, we obtain a mobility proportional to In^ [kpo). More 
accurately, one should use an expression for the mobility obtained by Andc®- (see Eq.(3.27) 
and Fig. 5 of that paper) but with the replacement of Z by Z* when calculating the strength 
of the Coulomb interaction. 



VI. SCATTERING OF DIRAC FERMIONS BY SHORT-RANGE IMPURITY PO- 
TENTIAL 

Here we will discuss quantum relativistic effects in the electrons scattering by a short- 
range potential. It appears that in the case of graphene the contribution of such defects to 
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the resistivity is essentially smaller than for the conventional nonrelativistic two-dimensional 
electron gas. We argue that this conclusion should help our understanding of remarkably 
high charge carrier mobilities observed in graphene^i^i^. 

Let us consider the case of a small concentration of point defects (to be specific, we 
will call them impurities) with the concentration Uimp and the angle-dependent scattering 
cross-section cr (0) . Then the defect contribution to the resistivity p reads^ 

2 1 

P 



r(M 



e^v],N{Ep)T{kp)' 

/rf0^^(l_cos<^) (43) 
J d<p 



2-K 



where N [Ep) = 2kp/'KfiVF is the density of states at the Fermi level (taking into account 
the double spin degeneracy and two valleys), r is the mean-free path. Note that the product 
VpN {Ep) is proportional to kp = ^Jwn in is the electron concentration) for both ultrarela- 
tivistic and nonrelativistic two-dimensional electron gas and thus any essential difference in 
their transport properties can be related only to the scattering cross- sect ion. 

To determine the scattering cross section one has to solve the two-dimensional Dirac 
equation which, for the case of massless particles and radially symmetric scattering potential 
V {r) takes the from (cf. Ref.-^ for three-dimensional case): 

^ + ^/i ('■) - - 1/ (r)l a (r) = 0, (44) 

where / = 0, ±1, ... is the angular-momentum quantum number, gi (r) e*'*^ and fi (r) e*('+^)''^ 
are components of the Dirac pseudospinor; to be specific we will consider the case of electrons 
E = hvpk > 0. 

Modifying a standard scattering theory^ for the two-dimensional case one should try the 
solutions of Eq. (jH]) outside the region of action of the potential in the form 

gi{r) = A[ji{kr)+UHi'\kr 
h (r) = I A {kr) + t.i/S {kr)\ , (45) 

where the terms proportional to Bessel (Hankel) functions describe incident (scattering) 
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waves, 

da {(f)) _ 2 
d(t> nk 



il(f> 

l=—oo 



(46) 



The Dirac equation for ultrarelativstic particles (14^ has as important symmetry with 
respect to replacement / < — > g, I < — > — / — 1 which means ti = Thus, Eq. (H6|l can be 

rewritten in the form 

da (0) 2 



irk 



(47) 



J^ti cos[(/ + 1/2), 

1=0 

The back scattering (0 = vr) is absent rigorously, as was discussed in the section 2. 

For the simplest case of the potential V (r) = Vq at r < R and V (r) = at r > R, using 
boundary conditions of continuity of the wave functions at r = R one finds 

^ Ji jqR) Ji+i jkR) - Ji jkR) Ji+i jqR) 
h\^^ {kR) Ji+, {qR) - Ji {qR) Hil\ (kR) 

where q = {E — Vq) /hvp. For the case of small energy, E <^ Vq, kR <^ 1 which is typical 
for graphene one has 



tdk)^^-^^^^{^) (49) 



TTi Ji+i{qR) fkR\^^^^ 
(Tif Ji{qR) \ Y) 

and thus the s-scattering (/ = 0) dominates. Substituting Eq. (H9!) into Eqs. fHTI) and (l43ll . 
one finds the estimation for the impurity contribution to the resistivity p ^ {h/Ae'^) nimpR"^- 
This means that scattering centers with the radius of potential R of order of interatomic 
distances and small concentration are irrelevant for the electron transport in graphene giving 
negligible contribution to the resistivity. 

If we would repeat the same calculations for nonrelativistic electron gas^ we would obtain, 
instead of Eq.dlHD 

^ jk/q) Ji [qR) Ji+i jkR) - Ji jkR) Ji+, jqR) 
hI'^ (kR) Ji+, (qR) - {k/q) J, (qR) H^^^ {kR) 

where k and q are, again, wavevectors outside and inside the potential region. In this case 
the s-scattering phase vanishes for small energies not linearly but only logarithmically and 
one obtains much larger resistivity 

^~A^n\n^kpR) ^ ^ 

(note that the same estimation holds for the case of massless Dirac electrons in a particular 
case of the resonant scattering when Jg {qR) = 0). 
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So, the basic conclusion is that whenever any potential works as a resonant scatterer in 
the case of massive particles, the scattering is non-resonant in the case of massless Dirac 
fermions in graphene, except some special values of parameters. 

Conclusions 

The examples considered demonstrate that graphene provides an unexpected bridge be- 
tween the condensed matter physics and quantum field theory. The impact of the ex- 
perimental discovery of this material on different areas of physics and industry can't be 
overestimated. First of all, graphene is the first example of truly two-dimensional crystals, 
in contrast with numerous quasi-two-dimensional crystals known before. This opens many 
interesting questions concerning thermodynamics, lattice dynamics and structural proper- 
ties of such systems which are, however, beyond the scope of this paper. Furthermore, 
single-layer graphene provides first experimental realization of a two-dimensional massless 
Dirac fermion system. The analogy with the quantum field theory proved crucial for under- 
standing of graphene unusual electronic properties, such as anomalous QHE, absence of the 
Anderson localization, inefficiency of scattering by point defects, etc. The bilayer graphene 
has a very unusual gapless parabolic spectrum providing an example of the system with 
electron wave equation different from both Dirac and Schrodinger ones. 
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